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Based on the model proposed by Hilgenfeldt at al. [Nature 398, 401 (1999)], we present here 
a comprehensive theory of thermal radiation in single-bubble sonoluminescence (SBSL). We first 
invoke the generalized Kirchhoff 's law to obtain the thermal emissivity from the absorption cross- 
section of a multilayered sphere (MLS). A sonoluminescing bubble, whose internal structure is 
determined from hydrodynamic simulations, is then modelled as a MLS and in turn the thermal 
radiation is evaluated. Numerical results obtained from simulations for argon bubbles show that 
our theory successfully captures the major features observed in SBSL experiments. 
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I. INTRODUCTION 

Single-bubble sonoluminescence (SBSL or simply SL), 
first discovered in 1989, is a phenomenon of periodic light 
emission by an oscillating gas bubble trapped in the pres- 
sure antinode of a standing ultrasound wave in water (or 
other fluids) (see |li,i^] for detailed reviews on SBSL). The 
oscillating bubble is stable enough to survive many days 
through billions of acoustic cycles while the produced 
flashes are highly regular and incandescent. The width of 
the emitted light pulse is around 10 — 100 ps with a peak 
power of the order 10 mW nil US- The light pulse has 
nearly a gaussian shape with a slight asymmetry, which 
is basically identical in the red and UV portions of the 
spectrum Q , and Hiller et al. further confirmed that 
the pulse width and the emission time were independent 
of wavelength. However, in an interesting twist, Moran 
et al. l6|| demonstrated that the pulse width did exhibit 
a mild dependence on wavelength at 3°C, but not at 
24 °C. Besides, the power spectrum of the emitted light 
was found to be broadband without any characteristic 
line, decreasing from the UV portion towards the red in 
a way that bore a resemblance to a blackbody spectrum 

SBSL has become an intriguing topic and an arena for 
experimentalists and theorists alike since its discovery. 
Numerous attempts have been made to study the bub- 
ble motion using classical bubble dynamics and sophis- 
ticated computational fluid mechanics (CFM) [gl llOlllll 
_13, 14, 15, 16] . It is generally believed that the bub- 
ble is heated to temperatures of tens of thousands and 
shock waves and plasma could be generated during the 
contraction of the bubble. Various proposals have been 
put forward to explain the light emission mechanism, in- 
cluding surface blackbody radiation H 0,0,0], neutral 
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and ion bremsstrahhm e; [T^ ITfll Hdi l2ll ] , collision- 

induced emission l23j]. quan tum vacuum radiation 
[2^ ] , conflned-electron mode l l25| , proton-tunnelling [2^ , 
and nuclear fusion I23, 12M [23 . Some of these propos- 

als, e.g. II S 111 El 111 El 111 E S III S 13 

attributed light emission in SBSL to the high tempera- 
ture attained in the bubble and are classifled as ther- 
mal radiation schemes in this paper. While qualita- 
tively reproducing the spectra detected in SBSL, most 
of such thermal schemes failed to explain why the pulse 
width is wavelength-independent as measured in some 
experiments 0, |1|. Owing to this major drawback of 
thermal schemes, researchers were forced to look into 
other non-thermal and more exotic models (see, e.g. 

mm mm El). 

To reconcile the success and the drawback of thermal 
radiation schemes, Hilgenfeldt, Grossmann and Lohse 
pol I21I ] took the finite opacity of the bubble into con- 
sideration and obtained a wavelength-independent pulse 
width for argon bubbles. The impact of Hilgenfeldt et 
aZ.'s work is huge and, to some extent, resurrects the 
thermal radiation scheme. However, it is worthwhile to 
note that the approach adopted in (23, i21i] is deemed a 
simplified version of the emission mechanism of SBSL as 
several p hysi cal processes have not been included in the 
proposal |30j. For example, the sonoluminescing bub- 
ble is modelled as a uniform one and the fluid dynam- 
ics inside the bubble has been neglected from the outset 
mm]. Besides, the Kirchhoff 's law used in [20II2H to 
evaluate thermal emissivity of the bubble has completely 
ignored the wave nature of light. 

In addition, existing hterature (e.g. H m, HI Ell) 
demonstrates there exists a gap in the theoretical treat- 
ment of light emission mechanism of sonoluminescence, 
in that blackbody radiation and thermal bremsstrahlung 
are often ascribed as separate possible mechanisms of 
SBSL. Often, Planck's formula for blackbody radia- 
tion 1^ Is^ and absorption coefficients for thermal 
bremsstrahlung in vacuum [20I \2l\ are applied separately 
in these cases, and the question of whether the bubble 
is opaque enough to demonstrate blackbody radiation 
is argued in a rather hand-waving manner by compar- 
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ing the photon mean path with an estimated size of the 
hght-emitting region. In our view, this is because of 
the lack of a single theory which can take account of 
both mechanisms in a finite-sized environment (i.e. the 
bubble) in a unifying manner. In a recent experimen- 
tal paper [s^ . Flannigan et al. demonstrated conclu- 
sively the existence of a plasma state inside the bub- 
ble, and hence thermal bremsstrahlung is an inevitable 
consequence because of the motion of the electrons and 
ions. We emphasize that blackbody radiation and ther- 
mal bremsstrahlung are nothing but one single emission 
mechanism manifested upon the degree of optical thick- 
ness of the bubble; and to this purpose, in this paper we 
have developed a coherent theory unifying both aspects 
in the context of SBSL. Within the framework of this the- 
ory, when the bubble becomes optically thick enough, the 
thermal bremsstrahlung manifest itself asymptotically as 
blackbody radiation. 

First, we will consider thermal emission in a realistic 
sonoluminescing bubble that is non-uniform in tempera- 
ture as well as density. To properly describe processes of 
thermal emission and absorption in a finite volume with 
a size comparable to the wavelength of light in a consis- 
tent manner, the generalized Kirchhoff's law is used in 
our paper |3J|. Second, state- of-the-art CFM is apphed 
here to evaluate the temperature and density distribu- 
tions in the bubble 1141 . (151] . Through such elaboration 
of the thermal radiation scheme |23, > '^^ succeed in 
obtaining an emission spectrum that agrees nicely with 
the experimental data as summarized in Ij and, in addi- 
tion, resolve the dilemma of whether the pulse width is 
dependent on the wavelength. 

The rest of this paper is organized as follows. We first 
present the generalized Kirchhoff's law in Sec. ^ and 
show in Sec . Illllth at it leads to the formula for emissiv- 
ity used in |20l l2ll | in certain limits. In Sees. IIVI and El 
we make use of the generalized Kirchhoff's law to derive 
the spectral radiance of a heated multilayered sphere. 
In Sees. IVII and IVIII respectively we furnish the plasma 
model and CFM used in the present paper. We then 
present relevant numerical results in Sec. IVIIII and con- 
clude our paper in Sec. IIXI 



We first review the statement of the 'classical' form 
of Kirchhoff's law (see, e.g. [s^)- Consider an isotropic 
absorbing medium with complex dielectric constant e = 
eR+iei, refractive index n = ^/e — nR+inj, and a typical 
size d. If the medium is in thermal equilibrium at tem- 
perature T and d is large compared with the wavelength 
of light A so that the geometric optics approximation 
holds, then the classical form of Kirchhoff's law relates 
the emission coefficient ri{w) and the absorption coeffi- 
cient k{uj) = 2niuj/c at (angular) frequency uj = 27r/ as 
El: 



where 



BUT) 



(2.1) 



(2.2) 



is the (frequency) spectral light intensity of a blackbody 
for one polarization, with c, h and ks being the speed 
of light in vacuum, the Planck constant h divided by 27r, 
and the Boltzmann constant, respectively. We remark 
that, besides the assumption A <C d, its application is 
justified only to a volume emitter which is optically thin. 

The generalized Kirchhoff's law (34] is a generalization 
of the classical form of Kirchhoff's law to all sizes d and 
optical thickness of a finite-size emitter, and can be de- 
rived from the Maxwell equations: 

V-D(r,t) ==p„(r,i), (2.3) 
V-B(r,<) = 0, (2.4) 

WxE{r,t) + ^^p^^O, (2.5) 
dD(r t) 

V X ff (r, t) ~ -^Q^ = Jn[r, t). (2.6) 

Here the electric field E, the magnetic induction B, the 
displacement field D and the magnetic field strength H 
are generated by the fluctuating charge density pn and 
current density J„, which are direct consequence of the 
random thermal motions of the charges and satisfy the 
fluctuation-dissipation theorem (see, e.g. [ssf): 



II. THERMAL RADIATION 

In SL, the high temperature reached inside the bub- 
ble ionizes the gas content, forming a partially ionized 
plasma [23, l2ll Issl Is^ . Besides having finite opti- 
cal thickness, the bubble has a size ~ 0.1 — 1 /ini near 
the instant of maximum compression, which is compara- 
ble to the wavelength of the emitted light ranging from 
200 nm to 800 nm. Hence, the assumption of geometric 
optics is invalid. To properly take account of the effects 
of finite absorption, wave reflection and diffraction^ we 
employ in this paper the generalized Kirchhoff's law |34j . 
which yields the 'classical' form of Kirchhoff's law as an 
asymptotic limit, to compute the power spectrum. 



(J„,j(r,u;)J* ,,(r',w)) 
4a; 

= —e{u:,T)e^ei{r,u:)5{r~r% 



(2.7) 



Here Jn.iif,^) is the i-component of the Fourier trans- 
form of Jn (r, t), 



.(.,T).-coth^ 



(2.8) 



is the average energy of a quantum harmonic oscillator 
at temperature T, and eo is permittivity constant of vac- 
uum. 

In the presence of the thermal fluctuations p„ and 
Jn, all fields E, D, H and B are thermally fluctuating 
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quantities. The set of stochastic equations, (|2.3|I - H2.6|I . 
together with the fluctuation-dissipation theorem, com- 
pletely determine the statistics of the electromagnetic 
field of such a system and the relevant theory is termed 
the thermal electromagnetic theory (TET) ^] . 

The generalized Kirchhoff 's law derived from TET [s^l 
states that the spectral radiance Pa, defined as the power 
emitted per unit wavelength interval and unit solid angle, 
is given by the following formula |34| : 



Pxin) = BxiT) CTabs(n). 



(2.9) 



Here P\{n) is the spectral radiance in the direction of 
n = r/|r|, (Tabs(n) is the absorption cross-sectional area 
of the emitter for an electromagnetic wave illuminating 
the emitter from the direction of n, and 



Bx{T) = 



(2.10) 



is the spectral light intensity in each polarization for an 
ideal blackbody. In 12.9|l it is understood that crabs is 
the sum of the absorption cross-sectional areas for inci- 
dent light waves with two perpendicular polarizations. It 
is evident that (Tabs is equal to the power Pa dissipated 
in the emitter for an incident plane wave carrying unit 
energy flux, which can be obtained from the following 
volume integral over the emitter: 

Pa = d^r eoei{r,Lo)\E{r,u:)\\ (2.11) 

It is worthwhile to note that E{r,uj) the electric field 
developed inside the emitter might be strongly enhanced 
at certain frequencies due to resonance effects and thus 
carries non-trivial frequency dependence [36|. 

Besides the spectral radiance, the integrated power of 
the emitted light pulse is another quantity measured in 
SBSL experiments and is simply the integral of Eq. (|2.9(l 
over the wavelength. 



P = 



dXPx 



(2.12) 



III. GEOMETRIC OPTICS MODEL 

We now apply the generalized Kirchhoff's law sum- 
marized above to consider light emission in the uniform 
bubble model (UBM) as proposed in 20, 21]. We will see 
that the formula for the spectral radiance used in poLl2]| 
is only an approximate one that is valid only under cer- 
tain restrictions. 

To use a simple situation to elucidate the generalized 
Kirchhoff's law, we first calculate the power emitted from 
a slab. Consider a weakly absorbing slab illuminated nor- 
mally by a plane-polarized plane wave with unit flux. Un- 
der the assumption that the slab (with thickness L and 




FIG. 1: The configuration of a uniform absorbing sphere. 



area A) has dimensions sufficiently large so that multi- 
ple internal reflections can be neglected and the inter- 
nal electric field can be represented by a single decaying, 
travelling wave only, the internal power loss can be found 
from Eq. ITTTJ: 



Pa 



^(l-p„f) Ail^e--^ 



KC 



(3.1) 



Here Pref is the fraction of power reflected from the slab 
surface. For small absorption, the imaginary part of the 
refractive index nj ~ e//2, so 



Pa - 2(1 - p,ef)BA(r)(l - e-'^^)A 



(3.2) 



where the factor of 2 properly takes care of the two 
possible polarizations. This result is consistent with 
that obtained from the standard radiative transfer the- 
ory l33j32|, upon which the spectral radiance obtained 
in 120,123 is founded. 

In the following we explicitly derive the formula for the 
spectral radiance that was used by Hilgenfeldt et al. in 
[20I [2H from Eq. H3.2|l . Fig. shows a uniform sphere 
of radius P, which is divided into multiple thin cylindri- 
cal shells. Each of these shells, indicated by the dashed 
lines, subtends an angle oi it ~ 29 at the center, and 
has a length L{0) = 2Rcos9 along the a;-direction and a 
differential cross-sectional area dA = 2ttR^ sin 9 cos 9 d9. 
With the assumptions that (i) p^d — 0, (ii) effects of 
refraction and diffraction at the spherical interface are 
negligible, (iii) internal reflection is ignorable, and (iv) 
each of these shells can be considered as a slab with area 
dA and thickness L{9), the radiance per unit solid angle 
follows directly from Eq. (|3.2() is given by: 



Pa = y" dA2Bx{T){l-e-^') 
= 2t: R^BxiT) 



,-2kR 



1 



- 1 



kR 



2k^R^ 



.(3.3) 



Multiplying this by the total solid angle Att straightfor- 
wardly yields the formula Eq. (17) in Ref. |23|. As seen 
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total absorption cross-section of the MLS. Eqs. H3.3() and 
(|4.1|) are the equations that we use for computing the 
power emitted from the SL bubble. In Sect. |Vl we will 
evaluate Cabsj from a wave optics perspective. 



WAVE OPTICS MODEL 




FIG. 2; An A''-layered spherical system. The inner A'^ — 1 
layers (the bubble) are absorbing, while the outermost layer 
(the surrounding fluid) is transparent. 



from the derivation here, this formula is only valid un- 
der the assumptions mentioned above, through which the 
wave nature of light has been completely neglected. Such 
an emission model will be referred to as the geometric op- 
tics model (GOM) in the following discussion. As we will 

the conditions 



show in our numerical results (Sec. IVIIip 
for the validity of GOM do not generally hold in a real- 
istic SL model. 



IV. MULTILAYERED SPHERE 



To evaluate Cabsj" in an MLS, we first determine the 
electromagnetic field inside the absorbing sphere using 
the transfer matrix formalism applicable to multilayered 
configuration (see, e.g. [s^). 

Consider a circularly polarized plane wave illuminating 
an absorbing MLS (see Fig.[2J). The MLS is composed of 
A^ — 1 spherical shells, and the refractive index, the inner 
and outer radii of the j-th shell (j = 1, 2, . . . , A^ — 1) are 
Jij, ttj-i and flj, respectively (ag = and ojv-i = R are 
assumed). The extended medium surrounding the MLS 
(i.e. the bubble) is considered as the A^-th layer and has 
a refractive index un- 

The electric and magnetic fields of the incident wave 
are given by Eiac{r) = {x + iy) exp{ikz) and Binc('") = 
—iE{r)/c. Hereafter we assume that the wave has pos- 
itive helicity. The electric and magnetic field inside the 
j-th layer of the MLS ( j = 1, 2, . . . , A''), can be respec- 
tively expressed in multipole expansion as: 



As discussed in Sect. lIIII Hilgenfeldt et al. [20|,|21| have 
assumed UBM as well as GOM in deriving the spectral 
radiance. This simplification serves as an illustration of 
the essential ingredients in SBSL. Yet its validity has to 
be verified, and for realistic calculations comparable to 
the experiment, it is necessary to take full account of the 
hydrodynamics inside the bubble. In light of this, in the 
following we will employ CFM developed by Ho et al. [s^ 
to simulate the hydrodynamics of the bubble. It is then 
obvious that UBM breaks down in such situation and, 
instead, we can model the inhomogeneous SL bubble as 
a multilayered sphere (MLS) with a piecewise-constant 
configuration of temperature and refractive index (see 
Fig. [2J|. In realistic calculations, the number of layers is 
so large that MLS is able to mimic the continuous dis- 
tribution obtained from CFM. On the other hand, MLS 
also includes UBM as a special case where there is only 
a single layer. 

Consider an MLS (i.e. the bubble) having a piecewise- 
constant temperature profile T = Tj, j ~ 1,2, . . . , A^ — 1. 
As the absorption cross-section simply becomes the sum 
of the contribution from each layer, it follows directly 
from Eq. that: 



1=1 



rijk 



BAr) = -^E*VM2I + 1) 



1=1 



V X [/j^Vg] 
rijk 



(5.1) 



(5.2) 



where the vector spherical harmonics Y^^yu^^jj = 
-irxVYim{e,4>)/^Jl{l + l), with Yim{0,(t)) being the or- 
dinary spherical harmonics ^| . 



f(M) 
■IN 



For the extended medium where j 
ji{nNx), with ji being the spherical Bessel function of 
the l-th order and x = kr, Eqs. H5.1|l - H5.2() reduce to 
the ordinary multipole expansion of a plane wave. For 
j < N, fj'^\r) {fj. = E, M) are governed by the radial 
wave equation: 



(4.1) 



dr^ 



2,2 1(1 + ^) 



r/fHr) = 0, (5.3) 



where (labs.j is the contribution of the j layer to the and, in addition, they satisfy the standard boundary con- 



5 



ditions on electromagnetic fields imposed at r = aj: 



_d 

dr 



p(M) 



(E) 



d 

dr 



rfft\{r) 



d r 



dr 



(M) 



d 

dr 



(5.4) 
(5.5) 
(5.6) 
(5.7) 



where h^j^^ is the Z-th order spherical Hankel function of 

the first kind. Substituting Eq. (|^ into Eqs. 

and recasting them into matrix form, we can show that: 



oip.) 



u 



(5.9) 



In general fj^\r) can be written as: 

f^f\r) = a\''^ji{n,kr) + I3\''h['\n,kr), (5.8) 



where T^^"* and T^^-* are respectively the TE-mode and 
TM-mode transfer matrix at r = aj, given explicitly by: 



-(E) 



(E) 



(E) luW, 



(E) 



[h\^\njxj),ji{nj+ixj)] 



(5.10) 



,(M) 



■ 2 2 



W^'^ljh\'^{n,x,),h\'>in,+^x,)] 



(5.11) 



r 



Here = kaj, and for convenience we define the gen- layer j = 1 to the layer j 



N 



1, i.e. T(^) 



eralized Wronskian Wj '^i j [/, g] for TE and TM modes 



as: 



get 



T 



(m) 



From Eq. H5.15|l we immediately 



fg' - f'g, 
fg' fg 



+( 



1 



where ' = d/dx. 



]_J_g_ 

^+1 ^ ^ 



In short, the transfer matrix T^''^ can be written as 



(5.12) 



(5.13) 



(5.14) 



and 



We can now solve for the field coefficients 

/Sj'^'' using the regularity condition at r = and radiation 
boundary conditions at r = 00, leading to the following 
relation: 



ai' = (^(z^))-!, (5.16) 
^(M) ^ c'^P)a[^^ =C^^^ /A^'^l (5.17) 

The field coefficients of each layer a^'^'' , /jj^"* is now readily 
obtained by iteratively applying Eq. 1)5. 9|) , from j = 1 to 
i — N — 1, using the boundary conditions Eqs. (|5.16|) - 

Now we proceed to calculate the absorption cross- 
section of the j-th spherical shell. To this end, we first 
evaluate: 



riNca-- 



Rejy" dQEj ■ (r x B*)| , 



(5.18) 



which is directly proportional to the energy flux J^j cross- 
ing the j-th interface at r = aj. With Eqs. H5.1|l . (|5.2|l . 
we find, after some algebraic manipulations. 



1 

Pn 







(5.15) 



Here A^'^', B(^), C^^), D^'^' (without the subscript j) are 
the elements of the total transfer matrix T*^^) from the 



fc2 



1 = 1 



,(5.19) 
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where x is evaluated at r = aj. Following directly from 
energy conservation, the absorption cross-section of the 
j-th spherical shell, Cabsj is given by the difference be- 
tween J-'j and Tj-i'. 

CTabsj = ^] - (5.20) 

Eqs. H5.19|l - H5.20|l in conjunction with Eqs. (|2.12l) and 

(14. 1|) , with no simplifying assumptions based on geomet- 
ric optics and optical thickness, are the major results of 
the wave optics model (WOM) introduced here, which 
will be used in this paper to calculate the light emission 
of a SL bubble. 



VI. PLASMA MODEL OF MULTILAYERED 
SPHERE 

A. Formation of Plasma 

In order to consider light emission processes in an MLS, 
we have to find the complex refractive index rij and hence 
the absorption coefficient Kj of each layer. In the presence 
of the high temperature developed during the collapse 
of a SL bubble, the gas inside the bubble is parti ally 
ionized at the instant of hght emission T^, T^, "^ff, '21', 
1^6]. Here we adopt a simple collision-dominated plasma 
model in which the collision frequency is a constant 
dependent on the concentration and temperature of the 
plasma |37| . The refractive index is then given by: 



where tOp = -y^ N^e^ /meQ is the plasma frequency, 

nb = -j T-j — t:^ , (6.2) 

is the contribution to the refractive index due to the back- 
ground neutral atoms as given by the Clausius-Mossotti 
equation (see, e.g. j^), and a is the atomic polarizabil- 
ity, with e and m being the charge and the mass of elec- 
tron, respectively. For convenience we have suppressed 
the j-dependence of the relevant physical quantities in 
the above equations. In the subsequent discussion we 
will use CFM developed in 35] to determine the number 
densities of electron and atom, and iVg, as a function 
of time. 



B. Collision Processes in plasma 

Around the instant of maximum compression, the tem- 
perature of a sonoluminescing bubble could reach tens 
of thousands of kelvins, ionizing the gas content inside. 
While Ref. j23| showed that the fraction of Ar+ ions 
amounts to less than 1% of the bubble content using a 



uniform bubble assumption; Ref. '35'| took account of full 
hydrodynamics and showed that this fraction can be as 
large as 30%. However, both papers point to the fact that 
the bubble content becomes a partially ionized plasma 
with Ar being the dominant species. Under this con- 
dition, bremsstrahlung is thought to play a major role 
in the light emission mechanism; in particular, electron- 
atom bremsstrahlung is expected to be the dominant pro- 
cess compared to other bremsstrahlung processes |20ll21j . 

On the other hand, as the numerical results obtained 
from CFM showed the degree of ionization may be much 
higher than those from UBM, we have developed here 
the WOM to take account of effects due to reflection, re- 
fraction and diffraction as well. In accordance with the 
approach of WOM, one has to consider the total effective 
collision frequency v, which is the sum of contributions 
from electron-ion collision, electron-atom collision and 
electron-ion recombination. According to the initial and 
final states of the electron, the first two are also known 
as free-free transitions, the latter one as free-bound tran- 
sition. Of course, the direct product of such collisions 
is the emission of photons and the corresponding mech- 
anisms are electron-ion bremsstrahlung, electron-atom 
bremsstrahlung as well as electron-ion recombination. 
Therefore, one can easily see the difference, as well as 
the relation between our approach and that proposed in 
[2CI _2^1J . In the following we provide the formulas of the 
collision frequencies in these processes. 



1. Electron-ion collision 

The simplest picture describing electron-ion collisions 
and electron-ion bremsstrahlung is to regard them as in- 
dividual binary events so that collective phenomena do 
not enter. Under such an assumption the differential 
emission cross-section du/, which measures the proba- 
bility of light emission due to the scattering of a unit in- 
cident electron flux from an ion with charge Ze, is given 
by the well-known Kramer's formula (39ij . It is custom- 
ary to include quantum mechanical corrections to this 
classical formula as the Gaunt factor |33|. For free- free 
transitions, including the free-free Gaunt factor gs{f,v) 
gives dfff as: 

where / is the light frequency and v is the speed of the in- 
cident electron. Within the range of optical frequencies, 
the free-free Gaunt factor is usually of the order 1. 

The collisions between charged particles is formally de- 
scribed by the scattering cross-section cr(6'), measuring 
the probability of scattering at an angle of a unit in- 
cident electron flux from an ion. However, the cross- 
section used in transport theory to predict the scattering 
frequency is Utr, related to the former through the rela- 
tion fTtr = cr(^')(l ~ COS0), where the over-bar indicates 
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averaging over the scattering angle 0. Using the relation 
between the differential emission cross-section dcjf and 
the transport cross-section trtr 



2, ,2 



do-/ 8 e V 



df Siireoc^hf''''' 
the transport cross-section is found to be 
4^2 / Ze^ 



O'tr 



V3 v47^eo"^^'^ 



(6.4) 



(6.5) 



The collision frequency, defined as v — Nvatr, where N 
is the number density of the background species {Ni for 
ions or Nq for atoms), is then 



47r2 N.Z^e^ 
^ (47reo)2TO2'i;3 



5ff(/, i^)- 



(6.6) 



2. Electron-ion Recombination 

In electron-ion recombination, an electron is first cap- 
tured by an ion, forming a bound state with energy levels 
labelled by quantum number n. A photon is released in 
such a process and the differential emission cross-section 
is given by Eq. H6.3|l with the free-free Gaunt factor re- 
placed by the free-bound Gaunt factor, gfb(n, /, w). In 
addition to its dependence on the photon frequency / 
and the velocity v of the incident electron, the free- 
bound Gaunt factor is also a function of n and approxi- 
mately of the order unity in optical frequencies. Accord- 
ingly, the transport cross-section and the electron-ion re- 
combination collision frequency are given respectively by 
Eqs. (|6.5() and (|6.6() with gfh{n,f,v) replacing gff{f,v). 



3. Electron- atom collision 

An electron moving near a neutral atom can also ex- 
perience a short-range Coulomb field, emitting radia- 
tion commonly known as electron-atom bremsstrahlung. 
There is no simple theory to determine the corresponding 
transport cross-section as we are aware of; and common 
practice is to determine it from experiment with differ- 
ent incident electron energies. With good accuracy in the 
relevant range of electron energies for an argon SL bub- 
ble, O'tr has a linear dependence on the electron energy 
Ee = mv^ jl and Ctr = Ct-cmv^ /2 -\- dtr, with the empiri- 
cal constants Ctr — 0.1 m^ and dfc — —0.6 x lO^^^m-^ 
|20l |4^ . The electron-atom collision frequency is there- 
fore 



Vce. = Nov{ctrmv'^ /2 -t- dtr)- 



C. Effective collision frequency 



(6.7) 



Assuming local thermodynamic equilibrium prevails in 
the plasma (as in the case for SL) and weak damping 



{i/ uj), we proceed to calculate the effective collision 
frequency defined by i/"^*^ = v = Nvatr{v), here 
indicates the Maxwellian average: 



47r 

Y \ 2TTkBT 



m 



X / dz;i;4e-™^'/2fe«T(...)_ 



(6.8) 



Accordingly, the effective collision frequencies for 
electron-ion bremsstrahlung, electron-ion recombination 
and electron-atom bremsstrahlung are obtained as (for 
clarity we drop the superscript 'eff'): 



'^ci(rc) ^ ^ ( "3" 



27r 



(6.9) 



8\/2 / kT \ ^ 

l^ea = -5-^0 (SctrfcsT+dtr). (6.10) 

Here the exponential factor e:xjp{—hf /ksT) is commonly 
referred to as the Gillie exponential factor i 5ff (/, T) 
and gfb(/, T) are the velocity-averaged free-free and free- 
bound Gaunt factor respectively (44] : 



Wif,T) = / di?.9ff(/,i?)e-^/'=-^(6.11) 

f^B-L Jhf 



gfh{f,T) = 2a;i V ^e^"5fb(n, /,«). 



(6.12) 



In H6.12|) . Xn = Eion/n^kBT, with E'jon being the first 
ionization energy of the atom, and n* is the lowest level 
for which En* < hf. 

Note that, assuming t^g, gfb ~ 1, the above equations 
can be reduced to 

Mf,T) - 1, (6.13) 

oo 

gEif^T) = 2x1 ^-e-". (6.14) 

The summation in Eq. (|6.14|) can be further simplified if 
the photon energy is small compared with the ionization 
energy, so that the energy level n* is high in comparison 
with the groimd state, as is the case for argon Since 
the density of the levels increases rapidly with increasing 
rt, the discrete levels higher than n* can be replaced by a 
continuous spectrum and the summation replaced by an 
integration, and Eq. (|6.14|l simplifies to 



1. 



(6.15) 



As the absorption coefficient k is related to the effective 
collision frequency as: 



«=(^^)2^, 
w c 



(6.16) 
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we can find the absorption coefRcients corresponding to 

Eqs. (ra- ionii : 



free-free Gaunt factor and exponential cut-off factor: 



'^ei(rc) 



27r 



3 3mkBT ('iTreo)^kBTc^m 



S'ff(fb)(/, T), 



(6.17) 



(6.18) 



The absorption coefficient Koi(rc) in 

Ref. !20| differs from 
ours in two ways, (i) They differ by a factor of kBT/hf. 
However, this difference is not very significant as kgT 
hf in SL; (ii) For electron-ion bremsstrahlung, the free- 
free Gaunt factor and the Gillie exponential correction 
were neglected in Ref. 20]. As we will discuss later (see 
Fig. d) , this could result in a factor of 3 difference in the 
calculated spectra. 

In the present paper, we explicitly take account of all 
Gaunt factors and the exponential correction, in par- 
ticular we adopt the fitting formula proposed by Itoh 
et al. 45„ .46.J to compute the average free-free Gaunt 
factor. For the free-bound Gaunt factor, since the pho- 
ton energies 1.5 — 6.2 eV, corresponding to the wave- 
length 200 — 800 nm, is small compared with the ion- 
ization energy (i?ion = 15.8 eV for argon), Eq. (|6.15|) 
is still a good approximation and we retain it for com- 
puting the average free-bound Gaunt factor. As a re- 
mark, following Eqs. 1)6. 9|l . (|6.13|) and (|6.15l) we have 
j/ic/fci — e'*-^/*^^"^ — 1, hence electron-ion bremsstrahlung 
is more dominant over recombination when hf <C fc^T. 
For SL, however, both processes are important since the 
thermal energy is typically 1.7 — 4.3 eV. 

For the purpose of comparison, in the following dis- 
cussion we will use two different sets of formulas, respec- 
tively denoted by PI and P2 models, to calculate the col- 
lision frequencies. The PI model employs the free-bound 
Gaunt factor, but ignore the free- free Gaunt factor and 
the Gillie exponential cut-off factor: 



2(^) 
^ 3 ' 



2( 



2tt. 



knT ^ 1 fcnT 



A-KtoksT 



hf 



(6.19) 



y^ksT^ikBT 
m hf 



3 ' iTTeoksT 
X9^if,T), (6.20) 

Z^ea = ^Na{'^)HiCtrkBT + dtr). (6.21) 

3 7rm 

Through (|6.16() . it is obvious that this set of formulas 
for the collision frequencies are consistent with the set of 
absorption coefficients used in [20l l2l| . 

By contrast, the P2 model, derived earlier in this sec- 
tion, readily takes account of free-bound Gaunt factor. 



2(?)^A^.( 



Ze2 



AneakBT 



xe -'--W(/,r), 



-hf/kBT—(^ 

,27r,3,_ Ze'^ 



2{—)'!Ni{- 

3 AneokBT m 

xe-'^f^'^^gE{f,n 

^N,{'^)H3ctrkBT + d,r). 

3 nm 



(6.22) 

(6.23) 
(6.24) 



Finally we state the basic assumptions underpinning 
the above formulas for plasma collision processes: (1) 
the plasma is 'cold' meaning that the electron thermal 
velocity is negligible with respect to the phase velocity 
of the wave, wth ^ "^ph; (2) the plasma is in the weak 
coupling regime, i.e. the ions are weakly interacting dur- 
ing their thermal motions. It is customary to indicate 
the degree of coupling by the dimensionless ion-coupling 
parameter (see, e.g. H^): 



2„2 



r = 



Z'e 



AneoRionkBT' 



(6.25) 



with i?ion — (47riVi/3)~3 being the mean inter-ionic dis- 
tance. If r ^ 1, the system is said to be weakly coupled; 
on the other hand if F > 1, the system is in the strong 
coupling regime '47| . We find F < 1 in a typical SL bub- 
ble (Sec. IVIIip . hence the assumption of weak coupling 
is at least approximately satisfied. 



VII. COMPUTATIONAL FLUID MECHANICS 

In this section, we summarize the CFM used in 
the present paper, which was developed by Gheng et 
al. jlj, [13 and later extended by Ho et al. |^ to include 
the ionization and recombination processes. The model 
couples the Rayleigh-Plesset (RP) equation for the bub- 
ble wall with the Navier-Stokes (NS) equations for the gas 
(including all the charged species resulting from ioniza- 
tions), while independently solving the energy equation 
for the surrounding water. The number densities of the 
charged species are computed from the reaction rates ap- 
proach. The effects of viscosity, surface tension, equation 
of state (EOS), compressibility and thermal conductivity 
of the ambient liquid are also taken into account. 



A. Bubble-wall dynamics 

To account for the effect of liquid compressibility, a 
more robust RP equation that gives the bubble radius R 
as a function of time t is used CM HM IHl : 



1 - M 
1 + M 



RR- 



3-M 
2(1 4- M) 



R'^H, 



Ps {t') , tnH, 



Pa 



1 + M' 



(7.1) 
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Here = R/ci, with ci the speed of sound in the sur- 
rounding Hquid, M = R/ci, t' = t + t^, po is the ambient 
hquid density, and Ps{t') = —Pa sin(wt') is the sonic driv- 
ing pressure with frequency lo and amphtude Pa- Also, 
the enthalpy Hi and the speed of sound of the liquid and 
ci are given by: 



Hi 



ci 



m 

Po pi ' 

dpi ' 



(7.2) 
(7.3) 



This modified RP equation includes terms to first order 
in the Mach number M of the bubble wall and allows for 
a variable q. 

Combining Eqs. ((OJ, lf73|l with the EOS of the am- 
bient liquid in the modified Tait form ||49.] . 



Pi + B 
Po+B 



(7.4) 



yields the explicit forms for Hi and ci: 



Hi 



ci 



B Pq + B 



Pi 



Po 



njPi + B) 
Pi 



(7.6) 



where B = 3049.13 bar and n ~ 7.15 are valid for water 
up to 10^ bar. 

Eqs. (|7.1I) . (|7.5|l and l|7.6|l must be supplemented by the 
boundary condition at the bubble wall, namely, that the 
pressure Pi (t) on the liquid side of the gas- liquid interface 
differs from the pressure P{R, t) on the gas side of the 
gas-liquid interface by the effects of surface tension and 
the normal component of viscous stresses , 



P{R,t)-T, 



rr I r—R 



m) 



jftiR 
R 



2a 



d{pE) 
dt 

2{pE 



or 

-P)v 



P)v 



(7.10) 



j,2 Qj, 



VTrr + 



k — 

dr 



Here r, p, P, T, Trr, k e and E — e + t;^/2 are the 
radial distance from the center of the bubble, gas den- 
sity, radial velocity, pressure, temperature, normal vis- 
cous stress, coefficient of thermal conductivity, the inter- 
nal energy and total energy per unit mass, respectively. 

If, due to ionizations and recombinations, there exists 
Ns species inside the bubble, then Ng — 1 mass conser- 
vation equations for these species must be supplemented 
with Eqs. H7.8|l - H7.10|l . In Ref. the maximum ion- 
ization level of the gas atom is taken to be 5, making a 
total of A^s = 7 species inside the bubble. This is more 
than adequate for the present temperature range; in fact. 
Ho et al. have shown that even the second ionization 
level can be safely ignored in practice. Note, since the ion 
densities change due to ionizations and recombinations, 
source terms must be added to the right hand side of the 
conservation equations. 

For convenience, let fj be the mass fraction of Ar^'^ 
(with j = 0,1,2,3,4,5) or electrons (with j — e), so 
(7.5) that X]j=e ~ ^- Therefore pfj represents the mass 
density of an individual species. The number density of 
an individual species is related to its mass fraction by 
Tij — pfj/rrij, where ruj is the mass of an atom (j — 0), 
or an ion with a charge j {j = 1 — 5), or an electron 
(j = e). The mass conservation equations of the species 
is then given by 



dipfj) , d 



dt 



dr 



{Pfjv) 



2pfjV 



(7.11) 



Here, the extra term {Ss)j is the source term for pfj 
which arises from ionization and recombination pro- 
cesses. It depends on the net rate of change of the number 
density of the species hj through: 



(7.12) 



For the ions {j = 0, 1, 2, 3, 4, 5), the net rate of change is 
given by: 



B. Hydrodynamics of gas 

The conservation of mass, momentum and energy for 
the gas flow in the spherical bubble is described by the 
compressible NS equations. They can be rewritten into 
a conservative form with source terms as: 



dp^d_ 
dt dr 



(pv) 



2pv 



(7.8) 



djpv) 
dt 

2pv^ 



d , 2 



P) 



r 



]_d_ 

2 Qj, 



■) + 



(7.9) 



(7.13) 



where aj™^^^ 



and 



are the rates of ion- 



ization, radiative recombination and three-body recom- 
bination of particles with a charge of j, respectively. The 
formulas for these rates can be found in |5(]l |. 

The net rate of change of the number density of elec- 
trons is simply given by charge conservation. Now that 



(7.14) 
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taking time derivative and multiplying both sides by p 
gives: 

5 

(5.)e = -5](5,),. (7.15) 

3=0 



C. Equation of State of gas 

The hydrodynamics of the bubble is certainly affected 
by the EOS. The most widely-used van der Waals EOS 
can be modified to take into account the ionization pro- 
cesses [Tgj i: 

p=(t- + -]r^^ (7.16) 

\ ^ rnj me J I - bp 
e = ^fcBT| + A] +A;sVV^T„ (7.17) 

2 I Too TOe / ^ ^ TOi 

\j=0 J " I J = l 2=j 

where Tj is the ionization energy of an ion with charge 
j — 1, and b the excluded volume. This EOS is denoted 

as MVEOS. 

The physical meanings of the MVEOS, Eqs. ifTT?^ and 
(|7.17|) . are manifest. The total pressure P is the sum of 
the contributions by different species, which are sepa- 
rately taken into account in proportion to their abun- 
dances. The internal energy e of the gas consists of both 
the thermal energy (the first term) and ionization energy 
(the second term). 



D. Energy transport in the liquid 

The changes in the liquid temperature T/ is accounted 
for with the assumption that the liquid compressibility 
and viscosity do not affect the heat transfer process be- 
tween the bubble and the surrounding water. As such, 
the energy equation for the water is: 

^+y^^ = Di-—(r^^] (7 18) 

dt dr r'^ dr \ dr J ' 

where vi and Di are the velocity and thermal diffusion 
coefficient of the liquid, respectively. 



VIII. NUMERICAL RESULTS 

In order to elucidate the significance of individual 
physical factors affecting SBSL, in the following wc 
present and compare numerical results obtained from 
simulations constructed with different models. Specifi- 
cally, we consider (i) UBM versus CFM model; (ii) COM 
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FIG. 3; A plot of the bubble radius R versus time t, which is 
normalized with respect to the acoustic period Td- 

versus WOM; and (iii) the two plasma models, PI versus 
P2. 

The rest of this section is organized as follows: first we 
study how the effects of various physical entities, includ- 
ing plasma, wave and temperature, can affect the emitted 
light pulse using the simple UBM model. Then we use 
the CFM model that is more realistic to mimic SBSL and 
compare relevant numerical results with those of UBM. 



A. Emission in Uniform Bubble Model 

UBM here refers to the model used in [20l l2l|. 

where the RP equation assumed incompressibility of the 
surrounding liquid and a variable polytropic exponent 
7(i?, R, T) was used to account for effects of thermal con- 
duction. However, instead of using the fitting formula in 
[20| for computing the polytropic exponent ^{R,R,T), 
we employed the formula proposed in [sj in the simu- 
lations. As in Ref. we studied the oscillations of a 
bubble with ambient radius Rq = 5.0 /im, subjected to 
an ultrasonic wave with / = 20 kHz and Pa = 1.3 atm. 
Figs. El El respectively show for the UBM the time evolu- 
tions of the bubble radius R over one acoustic cycle and, 
in the vicinity of the maximum bubble compression, the 
radius and the polytropic exponent 7. The number 
density of argon atoms and the temperature near the in- 
stant of minimum radius are shown in Fig. [SI Here we 
remark that the temperature profile is slightly different 
from that in [23| due to the difference in the formulas 
for 7(i?, R, T). The profile remains essentially the same, 
but the peak temperature in our result is about 3000 K 
higher. The temperature and density profiles are used as 
inputs in our calculations of the spectral radiance, from 
which other light emission properties (e.g. pulse shapes 
and FWHM) are obtained. 
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10' (t-U)/T, 

FIG. 4: The bubble radius R (upper panel) and the polytropic 
exponent 7 (lower panel) are plotted against the normalized 
time 10^ (t — tmin)/7d, where tmin is the moment at which R 
attains its minimum value of about 0.7/im. 
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FIG. 5: Time profiles of temperature T (upper panel) and 
number density of argon atom no (lower panel) near the in- 
stant of minimum bubble size. 
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FIG. 6: Free-free (solid line) and free-bound (dashed line) 
Gaunt factors multiplied by the Gillie exponential factor at 
the instant of minimum bubble size. 
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FIG. 7: (Color online) The spectral radiance Pa is plotted 
against the wavelength A for the following models: GOM-I-Pl 
(grey-dashed line), GOM-I-P2 (grey-solid line), WOM-fPl 
(dark-dashed line), and W0M-I-P2 (dark-solid line) 



1. Plasma and wave effects 

As mentioned previously, we employed two plasma col- 
lision models PI and P2 in the simulations. The differ- 
ences in these two models originate from the free-free 
(or free-bound) Gaunt factor and the Gillie exponential 
cut-off factor exp(— /i//fcsT). The exponential factor is 
usually close to unity in the Rayleigh-Jeans limit where 
hf <C fcsT, but for the case of SL, /cbT is of the or- 
der of a few eV and is within the range 1.5 — 6.2 eV of 
the observed light spectrum, and hence is not negligi- 
ble. In particular we find, when the bubble is at mini- 
mum size and the temperature and density of its contents 
are also at their maxima, the free-free Gaunt factor and 
the free-bound Gaunt factors, multiplied by the expo- 
nential cut-off, result in a correction factor of order 0.1 
(see Fig. IS)). Fig. [T] shows the computed power spectra 
obtained from simulations constructed respectively with 
G0M/W0M-I-P1/P2, clearly demonstrating that the P2 
model indeed leads to an decrease in the radiance. 



Furthermore, we can observe in Fig. that the power 
is overestimated when GOM is used rather than WOM, 
which can be readily explained as follows. When a plane 
wave is incident on a WOM bubble, part of it is reflected 
or scattered from the boundaries and the remaining part 
is absorbed as heat. In light of Kirchhoff's law then, 
less absorption implies less emission. In contrast, reflec- 
tion and diffraction are neglected in GOM, resulting in 
an overestimated absorption and hence emission. There- 
fore, to achieve realistic power calculations comparable 
with experimental results, the effects of both the wave 
nature of light and the Gaunt factor correction cannot 
be neglected, consistent with the point we made earlier. 

In subsequent discussions we will employ two specific 
light emission models: the model proposed in Ref. |23| 
(GOM-KPl) and our present model (W0M-HP2), and it 
should be understood that all power computations em- 
ploying GOM are done with PI while those employing 
WOM are done with P2. 
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FIG. 8: FWHM obtained from UBM+GOM versus wave- 
length A. Lines with circles, squares, crosses, diamonds and 
triangles respectively represent the cases with a scaling of 0.5; 
0.75; 1 (i.e. no scaling); 1.50; and 1.75. 



FIG. 10: FWHM versus wavelength A calculated from 
UBM+WOM. Cases with T = 20°C and T = 2.5 °C are in- 
dicated by the solid curve and the dashed curve, respectively. 
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FIG. 9: FWHM obtained from UBM-I-WOM versus wave- 
length A. Lines with circles, squares, crosses, diamonds and 
triangles respectively represent the cases with a scaling of 0.5; 
0.75; 1 (i.e. no scaling); 1.50; and 1.75. 



2. Temperature effects 

For our case studied here using the UBM model, the 
maximum temperature achieved at the instant of min- 
imum bubble size is about 23000 K (Fig. |SJl. To study 
the effect of the interior temperature on light emission, 
we scale by hand the original temperature profile for 
Pa — 1.300 atm and Rq — 5.0 /im by some chosen fac- 
tors, say 0.3,0.5,0.75,1.5 and 1.75, while keeping the 
densities of Ar neutrals constant. Accordingly, two of 
the three input parameters to the light emission model 
(temperature and ion number density) are changed and 
one (the atom number density) remains fixed. We employ 
both GOM and WOM to calculate the power, and study 
the spectral variation of the FWHM calculated within 
100 nm wavelength windows as shown in Figs. |H1 and |51 

It is remarkable that a consistently smaller spectral 
variation of the FWHM is obtained when the light emis- 
sion model WOM is used. The unsealed profile is about 



constant over the range of wavelength considered as in 
Ref. [23, which employed UBM and GOM, and the pro- 
file still remains remarkably constant when the temper- 
ature was scaled down by a factor of 0.5 and 0.75 (to 
a peak temperature of about 10000 — 20000 K). Scaling 
down by 0.3 (to a peak temperature of about 7000 K) 
produced zero power output since the temperature was 
much lower than that required for ionization. Neverthe- 
less, scaling up by 1.5 and 1.75 (to a peak temperature of 
about 35000 — 45000 K) produces dramatic variation of 
the FWHM. Hence, we find that the spectral uniformity 
of the FWHM holds only when the SL bubble tempera- 
ture is restricted within a rather small range of moderate 
values. However, we remark that the temperature is un- 
derestimated in the UBM since local temperature rises 
were not taken account of. With a more realistic hy- 
drodynamic modelling, Ref. found the temperature 
should be several 10^ K higher. Thus, we expect that a 
larger FWHM spectral variation with increased driving 
pressure will be an essential realistic feature of SL. 

In addition, we have studied the effect of ambient wa- 
ter temperature on SL light emission, using WOM and 
the values of {Pa, Rq) extracted from the phase diagram 
in Ref. p], at Tq = 20 °C and 2.5 °C and at a driving 
frequency of / = 26.5 kHz. Here we observe from Fig. 1101 



To = 20 °C 


Pa (atm) 


1.275 


1.283 


1.292 


1.300 


Ro (/im) 


2.9 


3.2 


3.6 


3.9 


To = 2.5 °C 


Pa (atm) 


1.320 


1.350 


1.375 


1.400 


Ro (/im) 


2.0 


3.5 


4.0 


4.5 



TABLE I: Driving pressure Pa and ambient radius Ro at 
To = 20 °C and To — 2.5 °C that give stable sonoluminesc- 
ing bubble for / — 26.5 kHz. 



the general trend of a larger FWHM towards the red end 
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of the spectrum as the pressure is tuned up, and this 
increase is further enhanced at a lower water tempera- 
ture. This is consistent with the experimental findings 
of Moran et al. ^ and is readily explained. At a lower 
water temperature, the bubble can be driven harder [s^ ] 
so that, given a certain value of ambient radius, a larger 
driving pressure can be applied while maintaining bubble 
stability. The bubble collapses more violently under the 
larger pressure and hence the temperature of the bub- 
ble interior achieves a higher value, resulting in a larger 
FWHM spectral variation. Therefore, in effect, both of 
our observations under increased driving pressure and 
lower water temperature can be explained in terms of 
the higher temperature reached inside the bubble. 



B. Computational Fluid Mechanics Model 

1 . Hydrodynamics 

Now we employ CFM developed by Ho et al. [s^ 
that includes the effects of ionizations and recombina- 
tions to calculate the power spectra and pulse profiles 
using WOM, and compare the results with those obtained 
from joint application of UBM and GOM. The set of con- 
ditions that we employ is extracted from Ref. [23| and 
shown in Table |Il| where Tq = 20 °C,/ = 20 kHz and 
the dissolved gas concentration is 0.20%. Figs. ^2 and 
respectively show the computed results of UBM and 
CFM for a case with Pa = 1.325 atm, Rq = 4.7 /im. It 
is found that the maximum temperature obtained with 
CFM can exceed 5 x 10** K while that in UBM is less than 
3 X 10^ K. Accordingly, the number of Ar+ ion in CFM 
is much greater than that in UBM. 



Pa (atm) 


1.275 


1.300 


1.325 


1.350 


Bjo (/xm) 


2.6 


4.0 


4.7 


5.4 



TABLE II: Driving pressure Pa and ambient radius Ro at 
To = 20 °C that give stable sonoluminescing bubble for / = 
20 kHz. 



2. Light emission 

To apply WOM to CFM, which produces an inhomo- 
geneous profile of bubble temperature and number densi- 
ties, we approximate the resultant inhomogeneous profile 
by a layered one and use within each layer j the average 
values of the temperature, the number densities of atom 
and electron there. The absorption coefficient n{j), colli- 
sion frequency ^{j) and refractive index can accord- 
ingly be computed using these averaged values. This ap- 
proximation scheme allows for the application of WOM 
to the resulting multilayered spherical configuration. 
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FIG. 11: Temperature T, number densities of argon neutrals 
n(Ar), ions n(Ar''') and electrons n(e~) shown as a function 
of time, where t = is the instant of minimum bubble radius, 
for a UBM bubble with Pa = 1.325 atm, and Ro = 4.7 pm. 




r (nm) r (fim) 

FIG. 12: (Color online) Snapshots near the instant of min- 
imum bubble radius {t = 0) for the same quantities as in 
Fig. 111! plotted against radial distance, for a CFM bubble 
with Pa = 1.325 atm, and Ro — 4.7/im. Grey lines indicate 
times before zero while dark lines indicate times after zero. 



As the degree of ionization in CFM result is much 
higher than that in UBM, the difference in the optical 
properties of these two models is obvious. As shown 
in Figs. [T^ and [TH cj/ujp is reduced by a factor of 2.5 
in CFM as compared to UBM. In particular, in CFM 
uj/ujp is close to unity near the UV end, indicating that 
plasma collective effects may be significant in the short- 
wavelength regime. Also, both dispersion and absorption 
are considerably stronger in the CFM case, exhibiting a 
larger variation in nji and a larger nj than the UBM 
case. 

In Sec. IVI Bl we have made the assumption that the 
plasma is so tenuous that the Coulomb energy is much 
smaller than the average thermal energy of individual 
particles and the plasma behaves like an ideal gas. The 
values of the ion-coupling parameter F at the instant of 
maximum bubble compression for the cases studied here 
(see Table |^ are shown in Fig. ^1 For UBM we note 
that F ~ 0.7 for all four cases. In the CFM bubble, in- 
terestingly there is a clear formation of two regions: an 
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FIG. 13: Shown on the left are the real part ur and imaginary 
part ni of the refractive index, on the right are the ratios uj/u)p 
and I'/ajp versus the wavelength A for a UBM bubble with 
Pa = 1.375 atm and Rq — 2.6 /xm at the instant of minimum 
radius. In the graph i^/up versus A, the dotted, short-dashed, 
long-dashed and full lines respectively show I'ei, Vrc, fea and 
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FIG. 14: Same as Fig. ESI for a CFM bubble with Pa = 
1.325 atm, and Rq = 4.7/im. The quantities shown here are 
those of the innermost layer which occupies the inner 7.5% of 
the bubble radius for Pa = 1.325 atm hence is representative 
of the hottest and densest region of the bubble. 
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FIG. 15: (Color online) The ion-coupling parameter F is plot- 
ted against r, normalized by the minimum radius i?miii at the 
instant of maximal bubble compression. The values obtained 
from UBM (CFM) are indicated by grey (dark) lines. 



FIG. 16: Spectral radiance Pa versus wavelength A obtained 
from GOM+UBM-I-Pl. The experimental spectrum is ob- 
tained from Barber et al. Q 



inner core which is moderately coupled T 0.7 and an 
outer shell which is weakly coupled F ~ 0. In addition 
the effect of increasing the driving pressure is seen to 
increase the size of this moderately-coupled inner core. 
Since the degree of coupling is moderate and not too 
strong, we expect that the formulas used for the absorp- 
tion coefficients and collision frequencies based on the 
tenuous plasma assumption should still apply. However, 
on the other hand, if the driving pressure increases while 
maintaining the stability of oscillation, it is likely that 
the plasma might become a non-ideal one. The physical 
property of such dense non-ideal plasma is rather com- 
plicated and is beyond the scope of the present paper 

In the following we contrast data obtained respectively 
from GOM-hUBM and WOM-fCFM and specifically con- 
sider three different physical quantities, namely the spec- 
trum, the pulse shape and the FWHM of light pulses. 

In Figs. El and El we show the computed spectra 
using GOM-HUBM-hPl and W0M-hCFM+P2, respec- 
tively. One clearly sees the improvement (Fig. 1171) of our 
refined model, namely W0M-I-CFM-I-P2, that the calcu- 
lated spectral shape is much closer to the experimental 
results [il than the GOM-t-UBM-fPl model. In particu- 
lar, major improvement is seen in the UV portion of the 
spectrum. 

The calculated pulse shapes are shown in Figs. 1181 and 
I19l respectivelv. evidently the pulse shapes produced from 
W0M-I-CFM-I-P2 are more consistent with experimental 
data 0, 0, where the long-time tail was not observed. 
By contrast, as shown Fig. El a long-time tail appears 
in GOM-fUBM-f PI and is an undesirable feature. 

Fig. 1201 shows the calculated FWHM plotted against 
the wavelength for both cases. While it is clear that 
the FWHM increases with driving pressure as found in 
experiments 0, 0, [s^ , we also see that the FWHM re- 
mains nearly a constant over 200 — 800 nm only at a 
low pressure Pa ^ 1.275 atm even in the more realistic 
W0M-|-CFM-fP2 model (c.f. Sec. lVIHA2| . We remark 
that the results obtained by Gompf et al. ^ showing sim- 
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FIG. 17: Spectral radiance Pa versus wavelength A ob- 
tained from WOM-I-CFM-I-P2. The best-fit for the experi- 
mental spectrum Q is indicated by the solid line, in which 
Pa = 1.276 atm, _Ro = 2.7 ^m. 
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FIG. 18: (Color online) Normalized power versus time ob- 
tained from GOM-l-UBM-fPl. The solid, dashed, and dotted 
lines are respectively the total normalized power, the normal- 
ized powers in the red and UV regions. 
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FIG. 19: (Color online) Normalized power versus time ob- 
tained from W0M-|-CFM-(-P2. The solid, dashed, and dotted 
lines are respectively the total normalized power, the normal- 
ized powers in the red and UV regions. 



FIG. 20: FWHM of light pulse versus wavelength A obtained 
at different driving pressures: circle Pa — 1.275 atm, square 
Pa = 1.300 atm, diamond Pa — 1.325 atm and triangle Pa = 
1.350 atm. Dashed lines with empty symbols and solid lines 
with filled symbols are calculated from GOM-|-UBM-|-Pl and 
WOM-hCFM-|-P2, respectively. 



ilar pulse widths for the red and UV pulse were obtained 
under a driving pressure of Pa = 1.200 atm; which is 
smaller than the lowest pressure Pa — 1.275 atm we used 
and is therefore expected to show a constant pulse width. 
In particular the results of Moran et al. ^ showed the 
FWHM increases with wavelength at a low ambient wa- 
ter temperature 3°C, which, as remarked previously, also 
resulted in a larger driving pressure. Thus we emphasize 
that in general the notion of wavelength-independence of 
the SL pulse width is only correct at low driving pres- 
sures; and at higher driving pressures spectral dispersion 
of the pulse width shows up, and this can be simply ex- 
plained within our model. Either effects of higher driv- 
ing pressure or lower water temperature boil down to 
the consequence of higher bubble temperature. As the 
bubble becomes hotter, both optical dispersion and ab- 
sorption become more significant and hence the bubble 
becomes more optically opaque, approaching a blackbody 
surface emitter. Red light is then radiated for a longer 
duration than the UV since, throughout one cycle, the 
bubble can stay at a lower-temperature state for a longer 
duration. As a consequence, the pulse width increases 
towards the red end of the spectrum. In other words, 
the emitted light becomes more spectrally dispersive be- 
cause of the increased absorption and dispersion in the 
plasma medium. In particular, the absorption is highest 
fFigs.[T^and [Tl|l at the red end of the spectrum, resulting 
in a flatter pulse shape there (larger FWHM). Thus, be- 
sides using a lower ambient temperature as in Ref . Q , if 
a strong enough pressure is applied (while still maintain- 
ing bubble stability) the spectral variation of the FWHM 
would be an observable consequence in experiment. 



16 



IX. CONCLUSION 

In summary, we have proposed in the present paper a 
robust theory for optical emission in SBSL that properly 
takes into account of the wave nature and propagation of 
light in the absorptive plasma formed inside a sonolumi- 
nescing bubble in a self-consistent way. In addition, our 
theory can be applied to bubbles with inhomogeneous 
density and temperature profiles. The validity of our 
scheme was examined for the case of SBSL with argon 
bubbles. By introducing and implementing appropriate 
Gaunt factors and exponential correction in the collision 
frequencies; as well as the effects of optical thickness, 
scattering, reflection and diffraction, our light emission 
model successfully explains the major features (includ- 
ing power spectrum, pulse shape and FWHM) observed 
in SBSL experiments. In addition, the computed power 
spectra and pulse shapes are shown to be in excellent 
agreement with experimental results. 

Besides, we have also shown that the experimentally 
observed spectral independence of the FWHM at Tq = 
20 °C is ascribable to the relatively small temperature 
range (about 10000 — 30000 K) achievable in a collaps- 
ing SL bubble. Above this range the bubble behaves in 
the way of a blackbody surface emitter and the spectral 
variation of the FWHM should be more notable. In fact, 



as the driving pressure goes up, the temperature reached 
inside the bubble rises. Also, if the ambient water tem- 
perature is lowered at a fixed ambient radius -Rq, the 
driving pressure allowable for bubble stability extends to 
a larger value. This provides a theoretical explanation 
for why Moran et al. 6] found a spectral variation of the 
pulse width at T = 3 °C. 

Notwithstanding the above-mentioned achievements, 
the model developed in the present paper is only one of 
the many steps towards a better understanding of SBSL, 
which is an extremely complex phenomenon resulting 
from the subtle interplay of hydrodynamics, chemical re- 
actions, plasma physics and optics as well. Much more 
challenging problems, e.g. SBSL with inert gases other 
than argon and inclusion of water vapor in the hydrody- 
namic code, are still ahead for us. They are surely our 
goal of endeavor in the future. 
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